REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

22-07-2010  Final  Report  22  July  2009  -  22-Jul-1 0 


4.  TITLE  AND  SUBTITLE 

Light  Pulse  Interaction  with  Narrow  Defects  in  Fiber  Bragg  Gratings 


6.  AUTHOR(S) 

Professor  Carlos  Martel 


5a.  CONTRACT  NUMBER 

FA8655-09-1  -31 06 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5d.  TASK  NUMBER 


5e.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Universidad  Politecnica  de  Madrid 
Plaza  Cardenal  Csineros  3 
Madrid  28040 
SPAIN 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

EOARD 

Unit  4515  BOX  14 
APO  AE  09421 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT  NUMBER(S) 

Grant  09-3106 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 


This  report  results  from  a  contract  tasking  Universidad  Politecnica  de  Madrid  as  follows:  -Description  of  work 
The  proposed  research  work  will  be  carried  out  in  2  steps: 

(I)  Formulation  of  the  problem  of  light  propagation  on  a  FBG  with  a  narrow  defect  in  terms  of  the  ID  Maxwell-Lorentz  equations.  Analysis  and 
numerical  computation  of  the  reflection/transmission  coefficients  for  the  CME  for  different  defect  profiles. 

(II)  Analysis  of  the  pulse  trapping/releasing  characteristics  of  the  defects:  study  of  the  internal  defect  modes  and  its  excitability  characteristics 
via  light  pulses. 


15.  SUBJECT  TERMS 

EOARD,  Fibre  Optics,  Solitons,  Nonlinear  Optics 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

UNCLAS 

b.  ABSTRACT 

UNCLAS 

c.  THIS  PAGE 

UNCLAS 

17.  LIMITATION  OF 
ABSTRACT 
UL 


18,  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
OF  PAGES  A.  GAVRI ELIDES 

64  19b.  TELEPHONE  NUMBER  (Include  area  code) 

+44  (0)1895  616205 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39-18 


Light  pulse  interaction  with  narrow  defects 
in  fiber  Bragg  gratings 

July  21,  2010 


Final  report  to  the  European  Office  of  Aerospace  Research  and  Development 
Research  grant:  FA8655-09-1-3106 


Carlos  Martel  (carlos.martel@upm. 

Dcpto.  Fundamentos  Matematicos 
E.T.S.I.  Aeronauticos 
Universidad  Politecnica  de  Madrid 
28040  Madrid,  SPAIN 


Disclosure  of  inventions:  I  certify  that  there  were  no  subject  inventions  to  declare 
during  the  performance  of  this  grant. 


Carlos  Martel  (carlos.martel@upm.es) 

Depto.  Fundamentos  Matematicos 
E.T.S.I.  Aeronauticos 
Universidad  Politecnica  de  Madrid 
28040  Madrid,  SPAIN 


Table  of  Contents 


List  of  Figures .  4 

Summary  .  5 

Introduction  .  6 

Methods,  Assumptions,  and  Procedures .  9 

Asymptotic  scaling  of  the  defects .  12 

Narrow  localized  defects  in  FBG .  16 

Results  and  Discussion .  25 

Computation  of  the  transmision/reflection  coeffcients .  25 

Nonlinear  defect  modes  .  29 

Defect  mode  stability .  34 

Pulse  trapping  simulations .  38 

Numerical  simulation  of  light  propagation  in  a  FBG  with  a  defect .  51 

Conclusions .  55 

References .  58 

List  of  Symbols,  Abbreviations,  and  Acronyms  .  60 


1 


List  of  Figures 


Figure  1  One  dimensional  fiber  with  a  perfectly  periodic  variation  of  the  refractive  index.  9 

Figure  2  Refractive  index  dependence  on  the  defect  profile  f(x) .  11 

Figure  3  Asymptotic  scaling  of  the  profiles  of  a(x)  in  eq.  (10):  a)  uniform,  b)  small 
periodic  structure  with  period  2tt//3,  c)  slow  defect,  slow  modulation  of  the 
small  periodic  structure,  d)  narrow  localized  defect  with  amplitude  of  the  order 
of  the  small  periodic  structure,  d)  narrow  localized  defect  with  amplitude  of 

order  1  centered  around  x  =  0 .  13 

Figure  4  One  dimensional  fiber  with  a  narrow  defect  in  the  periodic  grating .  17 

Figure  5  Sketch  of  the  dispersion  relation  (18),  with  the  resonant  wavenumber  k  =  1 

shown .  18 

Figure  6  Sketch  of  the  different  asymptotic  regions  in  the  x  —  t  plane:  (i)  1  and  2,  away 
from  the  defect,  where  the  NLCME  applies,  and  (ii)  the  narrow  shaded  band  of 

the  defect  in  the  vicinity  of  x  =  0 .  21 

Figure  7  Sketch  of  the  scattering  problems  that  give  the  jump  conditions  at  the  defect: 

(a)  outgoing  wavetrain  to  the  left  with  amplitude  1,  and  (b)  outgoing  wavetrain 

to  the  right  with  amplitude  1  .  22 

Figure  8  Top:  defect  profile  f(x)  =  De~x  ,  with  D  £  [0,1].  Middle:  modulus  of  the 
transmission  and  reflection  coefficients.  Bottom:  phases  of  the  transmission 

and  reflection  coefficients .  26 

Figure  9  Top:  defect  profile  f(x)  =  De~x  ,  with  D  £  [—0.85,0].  Middle:  modulus  of 
the  transmission  and  reflection  coefficients.  Bottom:  phases  of  the  transmission 

and  reflection  coefficients .  27 

Figure  10  Top:  defect  profile  f(x)  =  —D{x  +  1/3 )e~x  ,  with  D  £  [0, 1].  Middle:  modulus 
of  the  transmission  and  reflection  coefficients.  Bottom:  phases  of  the  transmis¬ 
sion  and  reflection  coefficients .  28 

Figure  11  Light  intensity  distribution  of  a  DM  localized  at  a  defect  at  X  =  0 .  29 

Figure  12  Standing  GS  profiles  for  9  =  7r/4,  tt/2,  and  7 r.  H(X)  =  R(X)e^x\  modulus 

R(X)  (top  plot)  and  phase  4>{X)  (bottom  plot) .  30 

Figure  13  Construction  of  a  DM  joining  two  pieces  of  a  standing  GS .  31 

Figure  14  Transmission,  ct,  and  reflection,  cr,  coefficient  size  for  the  defect  mode  families 

DM1  and  DM2 .  32 


2 


Figure  15  Sketch  of  the  profiles  of  the  defect  modes  in  families  DM1  and  DM2 .  33 

Figure  16  Stability  of  the  DM1  family  as  a  function  of  X\  and  ip.  Red  dots  indicate 
instability.  The  solution  is  also  shown  for  three  representative  cases  (marked 
with  black  dots).  Top  plots:  spatial  profiles  of  |^4+|  (blue  line)  and  |A~|  (red 
line)  at  the  end  of  the  simulation,  and  the  profile  of  the  DM  used  as  initial 
condition  (black  line).  Bottom  plots:  temporal  evolution  of  the  norms  1 1  \Af  2 1||.  36 

Figure  17  Stability  of  the  DM2  family  as  a  function  of  X\  and  ip.  Red  dots  indicate 
instability.  The  solution  is  also  shown  for  four  representative  cases  (marked  with 
black  dots).  Top  plots:  spatial  profiles  of  |A+|  (blue  line)  and  |A_|  (red  line) 
at  the  end  of  the  simulation,  and  the  profile  of  the  DM  used  as  initial  condition 

(black  line).  Bottom  plots:  temporal  evolution  of  the  norms  2||| .  37 

Figure  18  Left:  moving  Gap  Soliton  at  t  =  0,  A+{A~)  thick  (thin)  line,  with  parameters 
<7  =  |,  k  =  1,  c  =  \  and  9  =  2.  Right:  domain  of  existence  of  the  GS  in  the 

c  —  7  plane  (the  dot  corresponds  to  the  GS  shown) .  39 

Figure  19  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.2  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  i/|7L+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |A-|  in  red  at  T  =  0  (thin  lines)  and  T  =  350 

(thick  lines) .  41 

Figure  20  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.3  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  \/|A+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |v4“|  in  red  at  T  =  0  (thin  lines)  and  T  =  300 

(thick  lines) .  42 

Figure  21  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.4  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  i/|A+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |7l+|  in  blue  and  |A“|  in  red  at  T  =  0  (thin  lines)  and  T  =  250 

(thick  lines) .  43 

Figure  22  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.5  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  1 ^4+ 1 2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |v4“|  in  red  at  T  =  0  (thin  lines)  and  T  =  200 

(thick  lines) .  44 

Figure  23  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.6  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  \/\A+ 12  +  |)4”|2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |A“|  in  red  at  T  =  0  (thin  lines)  and  T  =  200 

(thick  lines) .  45 

Figure  24  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.2  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  \/|A+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |v4_|  in  red  at  T  =  0  (thin  lines)  and  T  =  350 
(thick  lines) .  46 


3 


Figure  25  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.3  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  ^/|yl+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |A“|  in  red  at  T  =  0  (thin  lines)  and  T  =  300 

(thick  lines) .  47 

Figure  26  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.4  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  ^/|A+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |v4“|  in  red  at  T  =  0  (thin  lines)  and  T  =  250 

(thick  lines) .  48 

Figure  27  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.5  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  ^/|4l+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |A-|  in  red  at  T  =  0  (thin  lines)  and  T  =  200 

(thick  lines) .  49 

Figure  28  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.6  sent  from  the  left.  Bottom: 

X  —  T  diagram  of  the  light  intensity,  I  =  \/|A+|2  +  \A~  |2,  during  collision. 

Top:  profiles  of  |A+|  in  blue  and  |v4“|  in  red  at  T  =  0  (thin  lines)  and  T  =  200 

(thick  lines) .  50 

Figure  29  Discretization  of  the  spatial  domain  [—L,L\ .  52 

Figure  30  Spatial  profiles  of  the  solution  at  T=10:  \A+\  (blue  line)  and  |A+|  (red  line). ...  54 
Figure  31  Modulus  of  the  difference  between  two  increasingly  refined  solutions  at  T=10. 

Blue  line  for  A+ and  red  for  A~ .  54 


4 


Summary 

This  is  the  final  report  of  the  research  grant  FA8655-09-1-3106  issued  by  the  Euro¬ 
pean  Office  of  Aerospace  Research  and  Development.  This  research  effort  is  devoted 
to  the  analysis  of  the  weakly  nonlinear  description  of  light  interaction  with  narrow 
defects  in  fiber  Bragg  gratings.  In  particular  it  includes:  (i)  the  derivation  of  the 
appropriate  envelope  equation  description  that  takes  into  account  the  effect  of  the 
grating  defect,  starting  from  the  full  ID  Maxwell- Lorentz  equations,  (ii)  the  anal¬ 
ysis  and  the  numerical  computation  of  the  resulting  reflection/transmission  coeffi¬ 
cients  induced  by  different  defect  profiles,  (iii)  the  derivation  of  a  family  of  nonlinear 
standing  defect  modes  that  keep  light  confined  to  the  vicinity  of  the  defect,  (iv)  the 
development  of  a  numerical  code  for  the  simulation  of  the  time  evolution  of  the  com¬ 
plete  problem  of  light  propagation  in  a  fiber  Bragg  grating  with  a  sharp  defect,  (v) 
the  numerical  analysis  of  the  stability  properties  of  the  nonlinear  defect  modes,  and 
(vi)  the  numerical  exploration  of  the  pulse  trapping  characteristics  for  some  typical 
defects. 
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Introduction 

Data  is  carried  as  light  pulses  in  optical  communication  systems.  A  very  high  reduc¬ 
tion  in  the  communication  speed  rate  comes  from  the  need  to  perform  the  processing 
of  the  signal  electronically.  It  is  therefore  very  interesting  to  try  to  find  new  materi¬ 
als  and  optical  structures  that  can  perform  all  signal  processing  optically,  bypassing 
the  need  for  the  optical/electronic  interface.  To  this  end,  we  have  seen  in  the  past 
years  an  increasing  research  activity  devoted  to  the  study  of  optical  nonlinearities 
and  light  interaction  with  nonuniformities,  in  an  effort  to  look  for  new  physical  ef¬ 
fects  that  can  be  applied  in  the  design  of  efficient  all-optical  switching  and  storing 
devices. 

One  of  the  most  promising  optical  configurations  are  fiber  Bragg  gratings  (FBG). 
FBG  are  microstructured  optical  fibers  that  present  a  spatially  periodic  variation  of 
the  refractive  index  along  the  length  of  the  fiber.  The  reflection  produced  by  the 
grating  induces  a  nontrivial  coupling  between  the  forward  and  backward  propagating 
light  beams  when  their  frequency  is  at  or  near  resonance.  This  effect,  combined  with 
the  Kerr  nonlinearity  of  the  fiber  that  appears  for  high  intensity  pulses,  results  in  the 
very  particular  light  propagation  characteristics  of  these  devices,  which  makes  them 
very  interesting  for  various  technological  applications  that  range  from  fiber  sensing 
to  all-optical  communication  (see  [1]  and  [2]  for  a  general  description  of  FBG). 

FBG  support  the  propagation  of  the  so-called  Gap  solitons  (GS),  which  are  a 
family  of  localized  nonlinear  light  pulses  that  propagate  over  the  periodic  grating. 
GS  are  not  true  solitons,  since  their  dynamics  is  not  completely  integrable,  but  they 
can  be  stable  and  propagate  without  distortion  along  the  FBG.  Standard  solitons 
in  a  bare  fiber  come  out  of  the  balance  of  nonlinearity  and  dispersion  [3],  and  this 
balance  can  take  lengths  up  to  several  kilometers  to  be  realized.  In  the  case  of  GS 
this  characteristic  formation  length  can  be  reduced  to  a  few  centimeters.  This  is 
because  GS  result  from  a  completely  different  balance  that  includes  propagation  at 
the  group  velocity,  grating  reflection,  and  nonlinearity.  Furthermore,  GS  have  the 
amazing  property  that  they  can  travel  at  any  speed  between  zero  and  the  speed  of 
light  of  the  bare  fiber.  This  possibility  to  have  “slow  light”  or  “trapped  light”  inside 
small  FBG  devices  has  made  them  extremely  promising  for  the  manufacturing  of 
all-optical  buffers  and  storing  devices,  and  have  fueled  numerous  studies  on  GS 
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dynamics  and  stability  in  FBG  in  the  past  15  years  (see,  e.g.,  [4,  5,  6,  7]  for  reviews, 
and  [8]  for  a  recent  experimental  realization). 

More  recently  a  new  strategy  for  trapping  light  inside  FBG  has  been  proposed 
that  consists  of  introducing  a  defect  or  inhomogeneity  in  the  periodic  grating.  The 
idea  is  to  send  a  moving  GS  that,  after  interacting  with  the  defect  in  the  grating, 
ends  up  trapped  in  the  form  of  a  stable  pulse  of  standing  light  localized  around  the 
defect. 

The  possibility  of  trapping  a  GS  at  a  grating  defect  depends  on  the  particular 
details  of  the  defect  and,  basically,  two  asymptotically  different  configurations  have 
been  analyzed:  slow  defects  and  narrow  defects.  Here  we  call  slow  defects  to  those 
that  consist  of  a  modulation  of  the  periodic  grating  on  a  scale  that  is  much  larger 
than  the  period  of  the  grating.  And,  on  the  other  hand,  narrow  defects  correspond 
to  distortions  of  the  grating  intensity  with  a  typical  length  scale  of  a  few  periods  of 
the  grating. 

For  slow  grating  defects,  the  problem  of  GS  trapping  and  releasing  has  been 
analyzed  in  several  papers  in  the  recent  years  [9,  10,  11,  12,  13,  14].  All  these 
analysis  are  theoretical  and  numerical,  and  have  been  performed  within  the  scope  of 
the  Nonlinear  Coupled  Mode  Equations  (NLCME)  formalism.  The  NLCME  are  a 
set  of  two  couped  envelope  equations  that  describe  the  weakly  nonlinear  dynamics  of 
light  propagation  in  a  FBG  of  small  intensity  (see,  e.g.,  [15]  and  references  therein). 
The  presence  of  the  slow  modulation  of  the  grating  gives  rise  to  a  spatial  dependence 
of  the  linear  reflection  coefficients  of  the  NLCME,  which  are  constant  for  a  perfectly 
uniform  periodic  grating.  The  analysis  of  the  characteristics  of  a  given  family  of 
defects  starts  with  the  computation  of  the  defect  modes,  which  are  standing  light 
pulses  localized  at  the  defect.  The  defect  modes  are  first  computed  in  the  linearized 
version  of  the  NLCME  for  small  light  intensity  and  then  extended  to  higher  intensity 
values  using  nonlinear  continuation  techniques.  Once  the  nonlinear  defect  modes  are 
computed,  one  has  to  tune  the  defect  parameters  in  order  to  find  the  stable  ones, 
which  are  the  only  valid  candidates  for  captured  light  pulses.  And  finally,  the  transfer 
of  energy  from  an  incoming  GS  to  a  stable  defect  mode  is  typically  analyzed  through 
the  numerical  simulation  of  the  dynamics  of  the  complete  process  using  the  variable 
coefficient  NLCME.  Following  these  steps,  GS  trapping  defects  have  been  located  in 
the  above  references  that  can  give  rise,  at  least  within  the  scope  of  the  theoretical 
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model,  to  captured  standing  light  pulses  inside  a  FBG. 

The  case  of  narrow  defects  is,  we  believe,  much  more  realistic  from  a  manufac¬ 
turing  point  of  view  than  the  previous  slow  modulation  defect,  but  has  not  been  so 
far  studied  in  that  much  detail.  The  analysis  that  can  be  found  in  the  literature  use 
the  NLCME  with  some  Dirac  delta  terms  added  to  somehow  take  into  account  the 
localized  effect  of  a  narrow  defect  (see,  e.g.,  [16,  17,  18]).  These  terms  are  added 
in  a  completely  heuristic  way,  just  to  try  to  qualitatively  reproduce  the  dynamics 
of  the  process,  and  they  have  no  connection  with  the  actual  shape  of  the  grating 
distortion. 

The  purpose  of  this  research  work  is  to  take  a  closer  look  into  the  interaction  of 
light  with  this  kind  of  narrow  defects  in  FBG,  and  it  is  organized  as  follows. 

In  the  first  part  (Methods,  Assumptions,  and  Procedures)  we  begin  with  a  discus¬ 
sion  of  the  correct  asymptotic  scaling  that  the  localized  defects  must  have  in  order 
to  appear  as  a  first  order  effect  in  the  NLCME  formalism.  We  then  show,  in  the 
next  section,  that,  from  the  point  of  view  of  the  slow  scale  description,  a  narrow 
defect  manifests  itself  as  an  internal  interface  inside  the  perfectly  periodic  FBG. 
The  resulting  transmission  and  reflection  coefficients  across  this  interface  depend 
on  the  particular  shape  of  the  defect,  and,  from  the  full  ID  Maxwell  equations  for 
light  propagation,  we  derive  the  scattering  problem  that  determines  the  values  of 
the  coefficients. 

The  second  part  (Results  and  Discussion)  presents,  in  the  first  section,  some  nu¬ 
merical  solutions  of  the  scattering  problem,  showing  the  transmission  and  reflection 
coefficients  induced  by  several  simple  defect  profiles.  In  the  following  section  we 
introduce  a  family  of  localized,  standing  nonlinear  defect  modes  that  keep  light  con¬ 
fined  to  the  vicinity  of  the  defects.  These  defect  modes  depend  on  the  particular 
values  of  the  coefficients  of  the  defect,  that  is,  on  the  particular  details  of  the  narrow 
defect  profile.  In  the  two  subsequent  sections  we  analyze  numerically  the  stability 
of  these  nonlinear  defect  modes,  and,  for  two  representative  stable  defect  modes, 
we  explore  the  possibility  of  light  pulse  trapping  (i.e.,  the  possibility  of  exciting  the 
standing  defect  modes  by  sending  moving  pulses  to  collide  with  them).  The  final 
section  of  this  second  part  is  dedicated  to  present  the  details  of  the  numerical  code 
used  for  the  computation  of  the  temporal  evolution  the  system.  And  finally,  this 
report  ends  with  a  summary  of  the  obtained  results  listed  in  the  Conclusions. 


Methods,  Assumptions,  and  Procedures 


For  the  case  of  a  fiber  with  a  perfectly  periodic  grating  and  a  cubic  Kerr  nonlinearity, 
we  will  use  as  a  model  for  light  propagation  the  one-dimensional  Maxwell’s  equations 
[1]  for  the  electromagnetic  fields  evolution  together  with  an  anharmonic  Lorentz 
oscillator  model  for  the  polarization  (see  e.g.  [9,  19]  and  references  therein), 


dB  _  dE 

dt  dx  ’ 
dD  _  dB 

^ 0  dt  dx  ’ 

D  =  e0E  +  P, 

d2P 

^p2~dt +  ^  _  2Ancos(27nr/As))P  -  7 P3  =  e0yP. 


(1) 

(2) 

(3) 

(4) 


In  the  system  above,  the  electric  field  P,  the  magnetic  field  P,  the  dielectric  displace¬ 
ment  P,  and  the  polarization  P  are  scalar  fields  that  depend  on  time  t  and  on  the 
spatial  variable  x  that  runs  along  the  fiber.  The  permeability  and  the  permittivity 
of  the  vacuum  are  denoted  by  ji0  and  e0,  respectively.  The  characteristic  frequency 
flp  accounts  for  the  non  instantaneous  polarization  response  of  the  media,  An  and 
A g  represent  the  strength  and  the  period  of  the  grating,  that  is,  the  strength  and 
the  period  of  the  spatial  periodic  variation  of  the  refractive  index  of  the  fiber  (An 
measures  the  size  of  the  nonuniformities  of  the  refraction  index  relative  to  its  mean 
value  n0,  see  Fig.  1),  y  is  the  linear  polarizability  of  the  medium  (n^  —  1  +  x)  and 
7  >  0  is  the  coefficient  of  the  nonlinear  Kerr  effect. 


n  0 


X 


Figure  1:  One  dimensional  fiber  with  a  perfectly  periodic  variation  of  the  refractive  index. 

In  order  to  simplify  subsequent  calculations  it  is  convenient  to  make  the  system 
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(l)-(4)  non  dimensional  using  the  following  rescalings: 


B  =  vW(eo7  )B,  D  =  (1/^)D,  E  =  (l/y/(^y))E,  P=(l/yfy)P, 

x  =  (Xg/n)x,  t  =  ( Xg/cTT)i , 


here  c2  =  l/(e0po)  is  the  vacuum  speed  of  light.  After  dropping  tildes,  the  nondi- 
mensional  Maxwell-Lorentz  equations  (MLE  hereafter)  can  be  written  as 


dB  _  dE 

dt  dx  ’ 
d  D  dB 

dt  dx  ’ 

D  =  E  +  P, 

d2P 

w” 2— -  +  (1  -  2Ancos(2x))P  -  P 3  =  (nl  -  1  )E, 
F  dt- 


(5) 

(6) 

(7) 

(8) 


where  the  dimensionless  finite  time  polarization  response  frequency  is  now  u2  = 
Wtc2  7 r2).  Notice  that  we  have  rescaled  the  spatial  variable  x  to  make  wavenum¬ 
ber  of  the  grating  equal  to  2  (see  eq.  (8)).  As  we  will  see  below,  the  reason  for  this 
choice  is  that  the  wavet rains  that  resonate  with  the  grating  and  develop  along  the 
fiber  will  then  have  wavenumber  1.  We  have  also  used  the  rescaling  to  absorb  the 
vacuum  properties  eo  and  po  (the  vacuum  speed  of  light  is  now  equal  to  1  in  the 
rescaled  variables)  and  the  nonlinear  coefficient  7. 

The  periodic  structure  of  the  fiber  appears  in  the  spatial  dependence  of  the  po¬ 
tential  of  the  anharmonic  Lorentz  oscillator  model  for  the  polarization  response  of 
the  media  given  by  eq.  (8).  In  order  to  take  into  account  the  possibility  of  a  defect 
in  the  grating,  we  change  it  to 


+  (1  +  f(z))P  ~p3  =  M  -  1  )£,  (9) 

where  f(x)  is  a  real  function  that  accounts  for  a  more  general  spatial  nonuniformity. 

Note  that  for  the  simplified  case  of  linear  light  propagation  with  instantaneous 
polarization  response  ( top  — >  00), 


•  if  there  is  no  spatial  structure,  f{x)  =  0,  then  we  obtain  from  eq.  (9)  P  = 
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(nl  —  1  )E,  which  once  inserted  into  eqs.  (5)-(7)  gives 

2d2E  _  d2E 
0  dt 2  dx2  ’ 

and  we  recover  the  standard  equations  for  linear  light  propagation  in  a  bare 
fiber  with  uniform  refractive  index  no- 

•  Also,  if  f(x)  =  — 2Ancos(2a;)  with  |An|  <C  1,  then  we  obtain 


(rig  +  2Ancos(2x)(nl 


!)  +  •■■) 


d2E 

W 


d2E 


that  corresponds  to  a  perfectly  periodic  spatial  modulation  of  the  refractive 
index  of  small  amplitude,  proportional  to  An  -C  1,  like  the  one  depicted  in 
Fig.  1. 


•  And  finally,  narrow  defects  can  be  taken  into  account  just  by  allowing  f{x) 
to  exhibit  a  more  general  profile  localized  in  a  spatial  region  of  characteristic 
size  5d  ~  1  (recall  that,  with  the  nondimensionalization  that  we  have  used,  the 
period  of  the  grating  is  equal  to  2).  The  resulting  effective  refractive  index  is 
then  given  by 


n 


2 


1  +  f(x)/n2 

1  +  f(x) 


where  a  positive  (negative)  value  f{x)  results  in  a  local  decrease  (increase)  of 
n2,  see  Fig.  2. 


-1  0 

Figure  2:  Refractive  index  dependence  on  the  defect  profile  f(x). 
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Asymptotic  scaling  of  the  defects 

Before  going  into  the  general  MLE  formulation  we  discuss  the  different  possibilities 
for  the  asymptotic  scalings  of  the  defects  and  how  they  translate  into  the  resulting 
envelope  equations  using  a  much  simpler  model  that  retains  all  qualitative  features 
of  resonant  wave  propagation  on  a  periodic  grating  with  defects. 

Our  toy  model  is  the  following  simple  linear  conservative  oscillator 

d2E 

where  we  allow  the  local  wavenumber  a(x)  to  depend  on  the  spatial  variable  x. 
We  are  going  to  consider  different  a(x)  distributions  (see  Fig.  3)  and  look  for  the 
corresponding  slow  scale  envelope  equations. 


Case  0  For  the  trivial  case  of  uniform  a(x)  =  «o  ~  1  (Fig.  3a),  the  exact  solution 
of  eq.  (10)  is  a  pure  harmonic  with  wavenumber  a0 

E  =  A0eiaoa:  +  c.c., 

where  A0  is  a  complex  constant,  and  c.c.  stands  for  the  complex  conjugate. 

Case  1  If  we  now  add  a  small  periodic  structure  with  wavenumber  (3:  a2(x)  = 
«o(l+£  cos  {(3x))  with  0<e<l  and  0  <  «o,  (3  ~  1  (Fig.  3b),  then  its  effect 
is  to  introduce  a  slow  modulation  of  the  amplitude  A0  that  evolves  in  a 
much  longer  spatial  scale  X  =  ex.  This  slow  modulation  can  can  be  easily 
captured  using  standard  asymptotic  perturbation  techniques  [20,  21].  The 
solution  is  expanded  as 

E  =  Eo(x,  X)  +  eE\(x ,  A )  +  . . . , 

which  once  inserted  into  eq.(10)  gives,  at  order  £°, 
n  2  e 

q^2  "F  aoEo  —  0  Eo  —  Ao(A")elQO‘l  +  c.c.. 

At  order  e1  we  obtain 

a  2  rp  d  A  • 

+  alEi  —  —  [2iao-^elQoa:  +  a20  cos(f3x)A0ela°x  +  c.c.]  =  RHS, 
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Figure  3:  Asymptotic  scaling  of  the  profiles  of  a(x)  in  eq.  (10):  a)  uniform,  b)  small  periodic 
structure  with  period  c)  slow  defect,  slow  modulation  of  the  small  periodic  structure,  d) 

narrow  localized  defect  with  amplitude  of  the  order  of  the  small  periodic  structure,  d)  narrow 
localized  defect  with  amplitude  of  order  1  centered  around  x  =  0. 

where  RHS  stands  for  right  hand  side  terms.  In  order  for  E±  to  be  bounded 
in  the  fast  scale  x,  the  following  solvability  condition  must  be  satisfied 

lint  -i-  [  RHS  e~iaox  dx  =  0, 

L— >oo  2 L  J 
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which  gives  the  evolution  of  the  amplitude  Aq  in  the  slow  scale  X : 


Case  2 


Case  3 


2i« 


dAn 


=  0 


dA0  ot{ 

2'a°lx  + 1 


dX 
Ao  =  0 


for  /3  7^  2a0,  and 
for  (3  =  2ao, 


where  we  have  take  into  account  that  cos (fix)  =  (e1/3x  +  e~1/3a;)/2,  and  that 
a0,  (3  >  0.  Note  that  only  in  the  resonant  case  /?  =  2a0  there  is  an  order 
one  effect  of  the  small  periodic  structure  in  the  amplitude  of  the  solution 
A0.  For  the  non- resonant  case  (3  ^  2 a0,  the  resulting  equation,  at  first 
order,  is  the  same  as  in  Case  0,  where  there  were  no  periodic  structure. 


In  the  resonant  case,  the  addition  of  a  slow  scale  modulation  to  the  am¬ 
plitude  of  the  small  periodic  structure,  a2(x )  =  q;q(1  +  sF(X)  cos(2ao^)) 
with  0  <  £  C  1,  0  <  a0  ~  1)  and  |F|  ~  1  and  (Fig.  3c),  just  produces  an 
evolution  equation  for  Aq  completely  similar  to  that  in  the  previous  case 
but  with  a  coefficient  that  depends  on  X: 

2i«°^  +  ^ F(X)A0  =  0 

This  case  corresponds  to  the  slow  defects  in  FBG,  which  are  described  by 
a  NLCME  but  with  variable  linear  reflection  coefficients  that  inherit  the 
slow  modulation  present  in  the  grating  [9,  10,  11,  12,  13,  14], 

A  narrow  distortion  of  the  periodic  resonant  structure  can  be  modeled  as 
a2(x)  =  ag(l  +  e(l  +  f{x ))  cos(2a0a;))  with  0<£<1  and  0  <  ao  ~  1, 
where  the  real  function  f{x)  is  zero  everywhere  except  in  the  narrow  defect 
region  of  length  5^  ~  1  where  \f(x)\  ~  1,  see  Fig.  3d.  The  corresponding 
evolution  equation  for  Ao  can  be  obtained  proceeding  as  in  Case  1.  The 
right  hand  side  terms  of  the  equation  for  the  correction  can  be  written 
as 

RHS  =  —  [2ia0^-^-e'a°x  +  tto((l  +  f{x)))  cos(2a0x) A0e^a°x  +  c.c.], 
dX 


14 


and  the  solvability  condition  now  gives 


2i«i 


dA0  al 

'ix  + 

1 


OZn 


lim  —  I  f(x)  cos(2a0x)(A0ela°x  +  c.c.)e  la°x  dx  =  0. 

ZIj 


'  —L 


Since  f(x)  is  zero  everywhere  except  in  the  narrow  defect  region,  the  in¬ 
tegral  above  is  bounded,  and  the  limit  as  L  — >  oo  is  zero.  The  resulting 
equation  for  the  slow  modulation  of  Aq  is  the  same  as  in  Case  1,  and  thus 
a  narrow  defect  in  a  small  periodic  structure  as  the  one  depicted  Fig.  3d 
produces  no  effect  at  all  at  first  order. 


Case  4  The  only  way  for  a  narrow  localized  defect  to  be  felt  at  leading  order  is 
to  have  a  much  higher  distortion  amplitude;  of  order  one,  see  Fig.  3e.  In 
this  case  we  can  write  a2(x )  =  chq(1  +  f(x)  +  £  cos(2a0a;))  with  0  <  £  C  1, 
0  <  a0  ~  1,  and  the  real  function  f{x)  is  again  zero  everywhere  except 
in  the  narrow  defect  region  of  length  5d  ~  1  centered  around  x  =  0  where 
1.  Away  from  the  localized  defect  the  structure  recovers  its 
periodicity  and  the  envelope  equation  obtained  in  Case  1  describes  the  slow 
modulation  of  A0.  So,  if  we  call  A\  and  A2  to  the  first  order  amplitudes 
of  the  solution  in  the  regions  at  the  left  and  right  of  the  defect,  we  have: 

1  A  2 

2ia0—^-  +  -^-Ai  =  0  for  X  <  0,  and  (11) 

dX  2 

1  A  2 

2ia0-7Tr  +  -77^2  =  0  for  X  >  0.  (12) 

dX  Z 

From  the  point  of  view  of  the  slow  scale  X  the  defect  is  seen  as  a  sharp 
interface  of  zero  thickness.  In  order  to  get  the  jump  conditions  across  this 
interface  we  have  to  study  the  internal  problem,  which  is  given  by  eq.  (10), 

^ +  ag(l +  /(*))£  =  0,  (13) 

where,  in  first  approximation,  we  have  neglected  the  effect  of  the  small 
grating  inside  the  defect.  This  equation  has  to  be  solved  with  the  boundary 
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condition 


E  — >  e'aox  +  c.c.  as  x  — >  — oo, 

and,  as  an  output,  one  gets  the  behavior  to  the  right  of  the  defect 
E  — >  Qelao:I  +  c.c.  as  a;  — >  +oo. 

Once  this  scattering  problem  is  solved  we  can  write  the  jump  condition 
across  the  defect 

Ao  =  CdAi  at  X  =  0.  (14) 

This  condition  together  with  eqs.  (11)  and  (12)  completes  the  description 
of  the  slow  scale  dynamics  of  eq.  (10)  for  a  small  resonant  structure  with 
a  narrow  defect.  Note  that,  for  the  computation  of  the  jump  condition  co¬ 
efficient  q,  the  detailed  knowledge  of  the  defect  is  essential,  and  the  small 
periodic  structure  is  not  relevant  since  it  is  of  higher  order  as  compared  to 
the  order  one  amplitude  of  the  defect.  It  is  also  worth  to  mention  that  the 
problem  inside  the  defect  given  by  eq.  (13)  appears  to  be  as  complicated 
as  the  original  problem  (10),  and  one  could  be  tempted  to  conclude  that 
no  simplification  comes  out  of  this  asymptotic  perturbation  analysis.  Well, 
this  is  obviously  not  true,  since  in  problem  (10)  we  have  a  short  spatial 
scale  x  ~  1  that  has  to  be  integrated  accurately  to  very  long  distances  of 
the  order  x  ~  1/e  1  (i.e.,  X  ~  1),  while  in  the  long  scale  description 

given  by  eqs.  (11),  (12)  and  (14)  there  is  only  one  spatial  scale  X  ~  1  and 
its  integration  is  numerically  much  cheaper. 


Narrow  localized  defects  in  FBG 


Now  that  the  required  asymptotic  scaling  of  narrow  localized  defect  has  been  made 
clear  in  the  previous  section  we  can  go  back  to  our  original  FBG  problem.  For  the 
sake  of  completeness  we  repeat  here  the  MLE  formulation  for  a  FBG  with  a  narrow 
defect  that  we  are  going  to  use: 


d2[E  +  P)  _  82E 

d t2  dx2  ’ 

d2P 

=  -wj(l  -  2Ajicos(2x)  +  f(x))P  +  -  1  )E  +  u)-p‘. 


(15) 

(16) 
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These  are  just  the  nondimensional  equations  (5)-(8)  after  elimination  of  the  dielectric 
displacement  D  and  the  magnetic  field  B.  We  are  interested  in  deriving  the  envelope 
equations  for  the  long  scale,  weakly  nonlinear  dynamics  of  this  system,  \E\,  \P\  <C  1, 
for  the  case  of  a  small  amplitude  periodic  grating,  An  <C  1,  and  when  there  is  a 
narrow  defect  of  extension  5d  ~  1  given  by  the  localized  function  f(x)  with  order 
one  amplitude  (see  Fig.  4). 


Figure  4:  One  dimensional  fiber  with  a  narrow  defect  in  the  periodic  grating. 


Region  away  from  the  defect 


To  the  left  and  right  of  the  defect  the  grating  is  perfectly  periodic,  and  the  appropri¬ 
ate  envelope  equations  in  these  regions  are  the  standard  NLCME.  We  briefly  sketch 
here  the  derivation  of  the  NLCME  (for  a  detailed  account  of  this  process  see,  e.g., 
[15,  22]). 

The  linear  propagation  characteristics  of  the  system  (15)-(16)  are  obtained  by 
looking  for  uniform  wavetrain  solutions  with  spatial  wavenumber  k  and  frequency 

E(x,  t)  If  Ek 
P(x,t )  J  |  Pk 

in  the  linearized  version  of  the  system  (15)-(16)  for  the  bare  fiber  (i.e. ,  in  the  absence 
of  grating  or  defects).  The  resulting  dispersion  relation  is  of  the  form  (see  Fig.  5): 


eikx+iukt+  cc  } 


(17) 


Uk  —  ± 


\j 02  +  ^>o)/2  ±  \JW +■  ^p«o)2/4  -  (18) 
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with  associated  eigenvectors, 


Note  that  there  are  two  very  different  behaviors  for  large  wavenumbers:  one  domi¬ 
nated  by  the  finite  time  polarization  response  of  the  medium,  u>k  —■ ►  ±up  as  k  — >  ±oo, 
and  a  second  one,  u>k  — ►  ±k  as  k  — >  ±cx),  which  corresponds  to  propagation  like  in 
the  vacuum  with  small  polarization  effects. 


Figure  5:  Sketch  of  the  dispersion  relation  (18),  with  the  resonant  wavenumber  k  =  1  shown. 


The  weakly  nonlinear  response  of  the  system  when  a  small  periodic  grating  is 
present  can  be  described  with  two  counterpropagating,  slowly  modulated  wa vet, rains 


E(x,  t ) 


P(x,t) 

which  have  wavenumbers  and  frequency 


=  V0(A+{x,t)eix+iuJt  +  A-(x,t)e~ix+i0Jt)  +  c.c.  +  . . . ,  (20) 


k  =  ±1  and  uk  =  u)  =  J\ (1  +  u%nl)/2  ±  ^/(l  +  ufoffi/A  -  c (21) 
and  resonate  with  the  grating  (recall  that  the  grating  wavenumber  is  2).  The  eigen- 
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vector  Vq  is  given  by 


Fo  = 


u 

1  —  u;2 


and  the  weakly  nonlinear  level  of  this  approach  requires  essentially  that 


•  •  •  <C  \AXX\  -C  \AX\  <C  \A±\  -Cl,  •  •  •  -C  \Af\  <C  |A±|  -C  1  and  An  -C  1,  (22) 


that  is,  small  amplitudes  that  depend  slowly  on  space  and  time,  and  small  grating 
strength. 

Once  we  take  the  ansatz  (20)  to  the  system  equations  eqs.  (15)-(16)  and  apply 
solvability  conditions  to  ensure  the  boundedness  of  the  solution  in  the  fast  spatial 
and  temporal  scales,  we  arrive  to  the  well  known  NLCME  for  the  evolution  of  the 
amplitudes  A± 


A+ =  vgA+ +  wAnA~ +  A+(Ul\A+\2 +  u2\A-\2)  + (23) 
At  =  — VgAx  +  wAnA+  +  A  (wi|A  |2  +  ■U2|A+|2)  +  . . . ,  (24) 


where  vg  is  the  group  velocity 


dui 


va  = 


dk 


k= 1 


“'(w2  -  w2) 
^-ojI  ' 


w  is  the  coefficient  of  the  effect  of  the  grating 

.  cu(l  —  a;2)  2 

W  =  1  2(uja  —  oj2)L°1p'> 


and  the  nonlinear  terms  are  given  by 

,3cu(l  —  cu2)3 


u  i  =  i- 


2  (^4-^2)  P 


uj  and  u2  =  i 


,3cu(l  —  u2)3  2 


ol 


(cu4  —  LUp)  P 


(25) 


(26) 


(27) 


ff  we  introduce  the  slow  scale  L  ~  1/An  1,  and  rescale  space,  time  and  the 
amplitudes 

x  =  LX,  t  =  ( L/vg)T ,  ~AA  =  y/vg/(L\u2\)A±,  (28) 
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we  arrive,  after  dropping  tildes,  to  the  scaled  standard  NLCME 


(29) 

(30) 


A%  -  A\  =  inA~  +  iA+(a\A+\2  +  \A~\2), 

Ar  +  Ax  =  ikA+  +  iA-(a\A~\2  +  |.4+|2), 

where,  in  order  to  make  the  nonlinear  and  grating  coefficients  positive,  we  have 
changed  to  the  complex  conjugates  of  the  amplitudes.  The  scaled  grating  strength 
k  =  AnL\w\/vg  ~  1  is  always  positive,  and  the  nonlinear  coefficient  a  =  \. 

The  NLCME  (29)- (30)  give  the  slow  time,  long  scale  evolution  of  light  propaga¬ 
tion  in  a  FBG,  and  represent  a  balance  of  the  transport  at  the  group  velocity,  the 
effect  of  the  grating,  and  the  cubic  Kerr  nonlinearity.  This  hyperbolic  system  and  its 
more  remarkable  solutions  have  been  widely  studied  theoretically  and  numerically, 
see  e.g.  [23,  6,  15,  24,  22]  and  references  therein.  As  sketched  in  Fig.  6,  the  NLCME 
are  valid  in  the  perfectly  periodic  grating  regions  1  and  2,  to  the  left  and  right  of 
the  defect,  respectively.  In  order  to  complete  the  formulation  of  the  problem  we 
have  to  study  the  narrow  defect  region  near  x  =  0  to  obtain  the  appropriate  jump 
conditions  that  would  prescribe  the  value  of  the  outcoming  envelopes,  A^  and  , 
in  terms  of  the  incoming,  ones  A j~  and  A£,  see  Fig.  6. 

Jump  conditions  at  a  narrow  defect 

In  the  narrow  region  near  the  defect  x  ~  1,  the  profile  \f(x)\  ~  1  in  the  MLE 
(15)-(16),  and  the  problem  can  be  simplified  to 

d2(E  +  P)  _  d2E 

dt2  dx2  ’ 

d2  P 

-w  =  -<4(i  +  f(x))p+<4(T%-i)E, 

where,  in  first  approximation,  we  have  neglected  the  small  nonlinearities  and  the 
small  grating  effects.  In  order  to  match  with  the  resonant  wavetrains  in  regions  1 
and  2  away  from  the  defect  (see  Fig.  6),  we  look  for  solutions  of  this  linear  system  of 
the  form  ( E,P )  =  (a(x),  b(x))elujt.  After  eliminating  b(x )  we  end  up  with  a  second 
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Figure  6:  Sketch  of  the  different  asymptotic  regions  in  the  x  —  t  plane:  (i)  1  and  2,  away  from  the 
defect,  where  the  NLCME  applies,  and  (ii)  the  narrow  shaded  band  of  the  defect  in  the  vicinity  of 


x  =  0. 


order  ordinary  differential  equation  for  a(x)  that  can  be  written  as 


d  a 
dx  ax 


0 

„  -4-(«o +/(a0) 

■co2lA - 

^-(i +/(*)) 

Up 


1 

0 


a 


(31) 


Note  that  the  detail  of  the  periodic  grating  is  only  present  in  the  equation  above  in 
the  frequency  u,  which  is  selected  through  the  wavenumbers  k  —  ±1  that  resonate 
with  the  grating  (see  eq.  (21)).  As  it  is  shown  in  Fig.  5,  we  choose  ui  on  the  lower 
branch  of  the  dispersion  relation.  This  lower  branch  is  the  one  that,  as  tup  — >  oo, 
gives  the  frequently  used  instantaneous  polarization  limit. 

Away  form  the  defect,  as  x  — >  ±oo,  f(x)  — >  0,  and  the  system  (31)  becomes  a 
constant  coefficient  problem  whose  solution  behaves  as 


a(x)  ~  a+e*x  +  a  e  lx. 


(32) 


The  propagation  characteristics  of  the  bare  fiber  are  thus  recovered,  and  the  coeffi¬ 
cients  a+  and  a-  have  to  match  with  the  left  and  right  going  envelopes  A+  and  A~ 
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of  the  regions  1  and  2  in  Fig.  6. 

To  evaluate  the  jump  conditions  at  the  defect  we  have  to  solve  the  two  scattering 
problems  depicted  in  Fig.  7,  which  correspond  to  the  determination  of  the  amplitude 
of  the  wavetrains  e1^  that  we  require  at  x  — >  ±oo  in  order  to  produce  amplitude  1 
for  elx  as  x  — >  — oo  (Fig.  7a),  and  amplitude  1  for  e-lx  asx^  +oo  (Fig.  7b). 


Figure  7:  Sketch  of  the  scattering  problems  that  give  the  jump  conditions  at  the  defect:  (a) 
outgoing  wavetrain  to  the  left  with  amplitude  1,  and  (b)  outgoing  wavetrain  to  the  right  with 
amplitude  1 

The  boundary  conditions  for  the  problem  in  Fig.  7a  are 

x  — >  —  oo,  a  — >  elx  +  cre~lx,  and 

x  — >  oo,  a  — >  ryelx, 

which  correspond  to  a  wavetrain  of  amplitude  1  outgoing  to  the  left,  and  to  no 
outgoing  wave  to  the  right.  In  order  to  extract  the  exponential  components  of 
the  solution,  we  multiply  the  solution  by  the  corresponding  left  eigenvectors  of  the 
constant  matrix  of  the  system  (31)  that  is  obtained  for  x  — »  ±oo  (see,  e.g.,  [25]). 
The  left  eigenvectors  are  given  by 

[±i,  1]  for  e^x, 
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and  the  above  boundary  conditions  can  be  transformed  into 


(33) 

(34) 


x  — >  —  oo,  ia  +  ax  — >  2ie1,:,  and 

a;  — >  oo,  ia  —  a^  — >  0, 

which  are  much  more  convenient  from  the  point  of  view  of  the  numerical  solution  of 
the  problem. 

Once  the  scattering  problem  (31)  with  the  boundary  conditions  (33)  and  (34)  is 
solved,  the  reflection  and  transmission  coefficients  can  be  obtained  from  the  limits 

x  — >  —  oo,  ia  —  ax  — >  2icre”l3;,  and  (35) 

x  — >  oo,  i a  +  ax  — >  2ictela;.  (36) 

The  problem  (31)  has  real  coefficients  so,  for  any  complex  solution  [a,ax\,  the 
complex  conjugate  [a,  dx]  is  another,  linearly  independent,  solution.  On  the  other 

hand,  the  trace  of  the  matrix  of  the  system  is  zero  (i.e.,  it  is  conservative)  and 

therefore  the  determinant  of  the  fundamental  matrix  composed  by  the  two  solutions 
[a,ax\  and  [ d,ax ]  has  to  be  constant  (see,  e.g.,  [26]): 

d 

—  \CLClx  ClClxj  0, 
ax 

which,  for  the  scattering  problem  in  Fig.  7a,  gives  the  following  relation 

|ct|2  + |cr|2  =  1.  (37) 

This  equation  just  states  the  fact  that  the  light  propagation  model  given  by  eqs. 
(15)  and  (16)  is  conservative,  and  thus  all  the  light  intensity  sent  into  the  defect  by 
the  incoming  beam  has  to  come  out  in  the  reflected  and  transmitted  beams. 

Moreover,  the  solution  of  the  scattering  problem  in  Fig.  7a  [a,  ax\  and  its  complex 
conjugate  [a,  ax\  can  be  also  linearly  combined  to  obtain  the  solution  of  the  scattering 
problem  in  Fig.  7b  [a',a^]: 

cj.  1 

[a  ,  a  ]  [a,  ax]  T  [a,  ax] . 

ct  Ct 

And  the  transmission  and  reflection  coefficients  for  this  second  scattering  problem 
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can  be  readily  obtained  without  having  to  solve  another  boundary  value  problem 

c't  =  Ct  and  c'r  =  -cr-r-  (38) 

ct 

Note  that  the  magnitude  of  the  transmitted  and  reflected  beams  is  the  same  from 
both  sides,  even  though  the  defect  is  not  symmetric. 

Finally,  taking  into  account  the  above  expressions  for  the  coefficients,  the  resulting 
jump  conditions  at  the  defect  for  the  slow  amplitudes  Af  and  Af  can  be  written  as 

At  =  ctAt  +  crA7,  and  An  =  ctA7  —  cr—At.  (39) 

'  '  “  Ct 

It  is  interesting  to  mention  that  for  the  particular  case  of  a  symmetric  defect, 


f(x)  =  f(~x), 


eq.  (31)  remains  invariant  under  the  change  x  — »  —x,  and  the  two  scattering  prob¬ 
lems  in  Fig.  7  are  just  the  same  problem  but  observed  form  opposite  directions.  The 
transmission  and  reflection  coefficients  of  both  problems  must  therefore  be  the  same, 


c't  =  ct  and  c'r  =  cr, 

and,  using  the  second  relation  in  (38),  we  can  conclude  that  for  a  symmetric  defect, 

\(t>ct  -0Cr|  =  (40) 

where  <f>Ct  and  (j)Cr  are  the  phases  of  ct  and  cr.  Or,  in  other  words,  for  a  symmetric 
deffect  the  reflection  coefficient  must  be  of  the  form  cr  =  iact,  with  «6l. 
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Results  and  Discussion 


Collecting  the  results  of  the  previous  section,  the  appropriate  envelope  equations  for 
the  description  of  the  slow  scale,  weakly  nonlinear  dynamics  of  light  propagation  in 
a  FBG  with  a  narrow  defect  at  x  —  0  (see  Fig.  6)  can  be  finally  written  as: 

A+t  -  A+x  =  inA^  +  h4+(cr|,4+|2  +  |Ab|2), 

A±T  +  A^x  =  inAf  +  iA^[{a\A^\2  +  |kL^|2),  for  X  <  0, 

A+t  -  A+x  =  ikA2  +  iA+(a\At\2  +  \A2\2),  (41) 

A2t  +  A2X  =  ikA2  +  i^42  (crl^2  I  "F \A+\  ),  for  X  >  0, 
together  with  the  jump  condition  (39): 

At  =  ctAt  +  crA7,  and  A2  =  ctA7  —  cr—At,  at  X  =  0, 

Ct 

where  k  is  proportional  to  the  depth  of  the  grating  and  a  =  \.  The  complex 
coefficients  ct  and  cr  are  determined  by  the  defect  internal  shape,  and  must  always 
verify  |ct|2  +  |cy|2  =  1. 

Computation  of  the  transmission/reflection  coefficients 

The  computation  of  the  reflection  and  transmission  coefficients,  ct  and  cr,  for  a  given 
defect  profile  f(x)  requires  to  solve  the  second  order  ODE  (31)  together  with  the 
the  boundary  conditions  (33)  and  (34). 

In  order  to  do  this,  we  first  move  the  boundary  conditions  to  a  finite  station 
x  =  i^oo,  where  the  defect  profile  has  sufficiently  died  away,  and  we  compute  two 
linearly  independent  solutions  corresponding  to  the  initial  conditions: 

[a,ax]  i  =  [1,0]  at  x  =  — Xoo,  and  [a,ax]  2  =  [0,1]  at  x  =  —  Xw 

The  general  solution  of  the  problem  can  be  written  as  a  linear  combination  of  these 
two  solutions 

[a,ax}  —  ccifa,  a^\\  +  0:2(0, 

where  the  complex  coefficients  aq  and  o?2  are  given  by  the  following  linear  system  of 
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equations 


i&l(  ^oo)  (  ^oo)  1^2  (  #oo)  Qjx2  (  ^oo) 

cq 

2ie 

i^'l(*^'C5o)  ^rrl(^oo)  ^2(^00)  ^x2  (*£00) 

0L2 

0 

which  results  from  the  application  of  the  boundary  conditions  (33)  and  (34).  Once 
the  solution  [a,  ax]  is  computed,  the  reflection  and  transmission  coefficients  are  read¬ 
ily  obtained  from  eqs.  (35)  and  (36): 

cy  =  (ia(-Xoo)  -  ax(-a:00))e_1Xoo/2i, 
ct  =  (ia^oo)  +  ax  (x00))e_lx°°/2i. 


We  have  computed  the  coefficients  for  the  following  three  representative  families 
of  localized  defects 

f(x)  =  De~x2  with  D  e  [0,1], 

f(x)  =  De~x2  with  D  e  [—0.85,  0], 

f{x)  =  -D(x  +  |)e-x2  with  D  e  [0, 1]. 


f(x)  “ 


hi  h\ 


fcthcr  « 


D  ' 


Figure  8:  Top:  defect  profile  /( x)  =  De  x  ,  with  D  £  [0, 1].  Middle:  modulus  of  the  transmission 
and  reflection  coefficients.  Bottom:  phases  of  the  transmission  and  reflection  coefficients. 
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Figure  9:  Top:  defect  profile  f(x)  =  De  ,  with  D  £  [—0.85,0].  Middle:  modulus  of  the  trans¬ 
mission  and  reflection  coefficients.  Bottom:  phases  of  the  transmission  and  reflection  coefficients. 


The  first  one  is  symmetric  and  positive,  the  second  is  also  symmetric  but  negative, 
and  the  last  one  in  not  symmetric  and  changes  sign  (see  the  top  plots  in  Figs.  8,  9 
and  10,  respectively).  The  remaining  parameters  of  problem  (31)  are  taken  to  be 

n2  =  2,  and  lo2  =  4, 

and  a i2  is  given  by  eq.  (21)  using  the  negative  sign,  see  Fig.  5.  The  numerical  inte¬ 
grations  have  been  performed  in  MATLAB  using  a  4th  order  Runge-Kutta  scheme 
with  error  below  10-6.  The  boundary  conditions  are  applied  at  the  finite  stations 
±Xao  with  Xoo  =  6,  and  we  have  checked  that  increasing  x^  produces  changes  in  the 
coefficients  that  are  way  below  the  selected  tolerance  of  10“6. 

The  modulus  and  phases  of  the  resulting  transmission  and  reflection  coefficients 
are  plotted  in  Figs.  8,  9  and  10.  It  is  interesting  to  mention  that  for  the  symmetric 
positive  defect  (Fig.  8)  the  transmission  coefficient,  ct,  is  only  slightly  reduced  when 
the  amplitude  of  the  defect  is  increased.  In  other  words,  the  defect  remains  almost 
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Figure  10:  Top:  defect  profile  f(x)  =  —D(x+  1/3 )e~x2 ,  with  D  G  [0,1].  Middle:  modulus  of 
the  transmission  and  reflection  coefficients.  Bottom:  phases  of  the  transmission  and  reflection 
coefficients. 

transparent  as  the  amplitude  of  the  defect  becomes  more  positive.  On  the  other 
hand,  when  the  defect  is  negative,  ct  decreases  much  faster  with  the  amplitude  of 
the  defect,  that  is,  the  defect  is  much  more  opaque  than  in  the  positive  case;  for  a 
negative  amplitude  D  =  —0.85,  ct  ~  .80  (see  Fig.  9),  while  for  a  similar  positive 
amplitude  D  =  0.85,  ct  ~  .995  (see  Fig.  8).  This  is  in  perfect  agreement  with  the 
results  plotted  in  Fig.  2,  which  showed  a  much  higher  sensibility  of  the  refractive 
index  for  negative  values  of  the  amplitude  of  the  defect  f(x).  From  the  bottom 
plots  of  Figs.  8  and  9  it  can  be  appreciated  that  the  difference  in  phase  between  ct 
and  cr  is  always  |,  as  it  is  expected  for  a  symmetric  defect  (see  eq.  (40)).  For  the 
third  case,  when  the  defect  has  both  positive  and  negative  parts,  the  situation  is 
somehow  in  between  the  previous  two:  ct  decays  slower  than  in  the  negative  case,  but 
much  faster  than  in  the  positive  case  since  the  negative  part  of  the  defect  has  higher 
amplitude  than  the  positive  one,  see  Fig.  10.  Note  also  that  in  this  non-symmetric 
defect  case  the  difference  in  phase  between  ct  and  cr  does  not  remain  fixed  to  |. 
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Nonlinear  defect  modes 


Defect  modes  (DM)  are  standing  light  pulses  that  remain  localized  around  the  defect. 
They  correspond  to  permanent  in  time,  light  intensity  distributions  that  are  zero 
everywhere  except  in  the  vicinity  of  the  defect,  see  Fig.  11.  In  this  section  we 
construct  this  type  of  solutions  of  the  complete  problem  (41)  by  joining  together 
two  parts  of  a  Gap  Soliton. 


/  =  l^+l2  +  \A~{2 


0 


X 


Figure  11:  Light  intensity  distribution  of  a  DM  localized  at  a  defect  at  X  =  0. 


The  Gap  Solitons  (GS)  are  a  well  known  family  of  exact  solitary  wave  solutions  of 
the  NLCME  (29)-(30)  that  propagate  without  distortion  at  any  velocity  c  in  the 
range  |c|  <  1,  that  is,  from  0  to  the  speed  of  light  of  the  fiber  without  grating  (see, 
e.g.,  [27,  23,  28,  29,  15]).  For  the  particular  case  of  zero  propagation  velocity,  c  =  0, 
the  standing  GS  can  be  written  in  the  form: 


A£s(X,T)  =  -^H(X)e~^T 

Ags(X,T )  =  V^H(X)e-^r 


where 


7  =  cos  6 ,  with  0  <  9  <  7r, 

and  the  function  H  is  given  by 

H(x)  =  s[n9 _ 

v  1  a  cosh[A"  sill  0  +  i#/2] 

The  standing  GS  are  pulse  solutions  with  steady  modulus  and  temporal  frequency 


(42) 

(43) 

(44) 

(45) 
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K'y.  The  complex  function  H  that  gives  the  prohle  of  the  standing  GS  verifies 


H{X)  =  H(-X).  (46) 

It  has  maximum  modulus  at  X  =  0  and  then  decays  to  zero  monotonically  and 
exponentially  fast  as  A"  — »  ±oo.  In  Fig  12  we  plot  the  profiles  of  the  standing  GS, 
which,  for  fixed  k  and  er,  depend  only  on  the  parameter  9. 


Figure  12:  Standing  GS  profiles  for  9  =  7r/4,  7t/2,  and  n.  H(X )  =  R(X)e^x\  modulus  R(X) 
(top  plot)  and  phase  <j>(X)  (bottom  plot). 


As  sketched  in  Fig  13,  we  can  join  two  pieces  of  a  standing  GS  to  form  a  DM.  The 
GS  are  exact  solutions  of  the  NLCME  and  therefore  the  resulting  DM  verifies  the 
differential  equations  for  Af  and  Af  in  the  system  (41).  It  must  verify  also  the 
jump  conditions  (39): 

Af  —  Ct(Af  +  ao4j“),  and  A^  =  Q(Af  —  aAf),  (47) 

with  a  =  Cr/ct.  The  values  of  Af  and  Af  are  given  by  the  expression  of  the  standing 
GS  (42)-(43)  at  X  =  Xi  and  X  =  X2  (see  Fig  13),  and  once  inserted  into  the  jump 
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conditions  produce: 


-R^1  =  ct(-R2e^2  +  aR le"^1), 
i?2e‘i<fe  =  CtiRie-^1  +  aR2e ^2), 


where  H(X\)  =  i^e1^1  and  H(X2)  =  R2e1(^2.  From  the  equations  above  we  can 
easily  obtain  the  transmission  and  reflection  coefficients 


ct  =  cos  "0 


2  R±R: 


Cr  =  (cOS  t/? 


Ri  +  Ri 

TD  2  d2 
rt2  -Tti 


pi((<^i-02)-V') 


f2?  +  i$ 


—  isin'0)e^2^>1 


(48) 

(49) 


which  depend  on  the  limit  points  of  the  GS  pieces,  Ad  and  (41)Ad,  and  on  the  free 
parameter  —  |  <  -0  <  |,  and  verify  the  conservation  condition  (37):  \ct\2  +  |cr|2  =  1. 


Xi  x  x2  x 


Figure  13:  Construction  of  a  DM  joining  two  pieces  of  a  standing  GS. 


Generally  speaking,  one  can  use  the  expressions  above  to  obtain,  for  given  values  of 
Ct  and  cr,  the  corresponding  DM,  which  are  defined  by  the  parameters  Xi,  X2  and 
ip.  Note  that,  because  of  the  conservation  condition  (37),  prescribing  the  values  of 
the  complex  numbers  ct  and  cr  corresponds  to  imposing  only  three  real  conditions 
and,  therefore,  we  have  the  appropriate  number  of  equations  for  the  determination 
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of  the  three  real  parameters  of  the  DM. 

Instead  of  analyzing  the  DM  in  the  general  form  given  by  (48)-(49),  we  will  restrict 
our  attention  to  two  particular  families  of  DM  that  share  the  property  R\  =  i?2 
(i.e.  no  jump  of  the  modulus  across  the  defect).  The  transmission  and  reflection 
coefficients  for  these  two  families  take  the  simplified  form: 

•  DM1:  X\  =  X2l  R\  =  Ro,  and  0i  =  0 2, 

Ct  =  cos^e-1^’, 
cr  =  —  isin^e^2^1-^ 

•  DM2:  Xi  =  —  X2,  Ri  =  R2,  and  0i  =  —02, 

Ct  =  cos^e^2^1-^,  (52) 

cr  =  -isin0ei(2^“^.  (53) 

The  expressions  above  depend  only  on  two  parameters,  —  00  <  X\  <  00  and  —  f  < 
0  <  jj,  and  they  span  all  possible  cases  between  total  transmission  cr  =  0  and  total 
reflection  ct  =  0,  see  Fig  14. 


(50) 

(51) 


Figure  14:  Transmission,  Ct,  and  reflection,  cr,  coefficient  size  for  the  defect  mode  families  DM1 
and  DM2. 


The  profiles  of  the  defect  modes  in  family  DM1  are  sketched  in  Fig  15.  They 
correspond  to  non-symmetric  defects  and,  for  0  =  0,  they  reproduce  the  trivial 
situation  of  no  defect  at  all  (ct  =  1  and  cr=0).  If  we  apply  a  spatial  reflection 
symmetry  to  the  complete  system  (41):  X  — >  —X  and  Af  A^,  then  the  equations 
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remain  invariant,  ct  does  not  change  either,  but  cr  becomes  c'r  =  -cr[|.  For  the  DM1 
family  this  implies  that  the  problem  for  the  DM  for  a  given  Xi  is  the  same  as  for 
—Xi  and,  therefore,  it  it  enough  to  study  the  case  Xi  >  0. 


DM1 


DM2 


0  X 


Figure  15:  Sketch  of  the  profiles  of  the  defect  inodes  in  families  DM1  and  DM2. 


The  second  family  DM2  corresponds  to  symmetric  defects,  i.e.,  cr  =  -crf-,  and  the  re¬ 
flection  symmetry  does  not  provide  now  any  reduction  of  the  parameter  space.  This 
is  precisely  the  family  of  symmetric  DM  derived  in  [17]  using  a  localized  defect  for¬ 
mulation  based  in  the  addition  of  Dirac  delta  functions  to  the  standard  NLCME.  In 
[17]  the  dynamics  of  the  defect  modes  was  analyzed  including  only  a  local  disruption 
of  the  Bragg  grating  reflection  and  without  local  perturbation  of  the  refractive  in¬ 
dex.  This  situation  corresponds  in  our  formulation  to  the  particular  case  of  ij)  =  20i, 
which  gives  real  ct  and  purely  imaginary  cr. 
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Defect  mode  stability 

We  analyze  the  stability  properties  of  the  defect  modes  obtained  in  the  previous 
section.  We  fix  6  =  7r/4  so  the  standing  GS  used  to  generate  the  DM  is  stable  (see 
[28,  29]),  and  we  study  the  stability  of  the  complete  families  DM1  and  DM2  for 
— oo  <  X1  <  oo  and  —  |  <  i/j  <  |. 

For  each  value  of  X\  and  ip  we  compute  the  corresponding  coefficients  q  and 
cr,  and  the  associated  DM.  We  then  run  a  numerical  simulation  of  equations  (41) 
using  as  initial  condition  the  DM  with  a  small  superimposed  random  perturbation. 
If  the  DM  persists  in  time  we  consider  it  stable,  and  unstable  if  the  system  evolves 
away  from  it.  The  details  of  the  numerical  method  used  to  integrate  equations  (41) 
can  be  found  in  the  final  section  of  this  chapter,  and  the  particular  values  of  the 
numerical  parameters  are  as  follows:  spatial  domain  X  G  [—20,  20],  final  integration 
time  T  =  200,  size  of  the  perturbation  added  10"3,  time  step  AT  =  .01,  and  grid 
spacing  A  A"  =  .01. 

The  stability  results  for  the  family  DM1  are  presented  in  Fig  16,  where  the  red 
dots  in  the  plane  X\  —  -0  indicate  the  unstable  DM.  The  DM1  family  is  X\  <->  —X{ 
symmetric  and  so  is  the  resulting  instability  region  (up  to  the  approximation  of  the 
method).  Note  also  that,  as  |Ab|  grows,  the  defect  is  located  further  in  the  tail  of 
the  GS  and  becomes  less  relevant,  and  thus,  for  large  |Ab|,  we  recover  the  stable 
dynamics  of  the  GS  without  defect.  A  typical  stable  DM  of  this  family  is  represented 
in  the  sub-figure  1.  In  the  top  plot  we  show  the  initial  (T  =  0)  and  final  (T  =  200) 
profiles  of  the  solution,  which  are  basically  indistinguishable.  The  bottom  plot  of 
the  sub-figure  shows  a  log  plot  of  the  temporal  evolution  of  the  quantity: 

||M,i 

II  QT  II 

for  the  four  amplitudes  Af  and  Af.  The  modulus  of  the  DM  is  steady  and  therefore 
this  magnitude  indicates  easily  the  departure  from  the  DM  solution.  For  this  stable 
case  this  magnitudes  decay  from  the  initial  perturbation,  and  remain  always  very 
small  ~  10-6.  Two  regions  of  instability  can  be  appreciated  in  Fig  16.  In  sub¬ 
figure  2  we  plot  a  solution  inside  the  instability  region  near  ^  =  7t/2.  The  regular 
oscillations  in  the  bottom  plot  indicate  that  it  corresponds  to  a  Hopf  destabilization 
with  a  clear  temporal  frequency  of  about  T  ~  8.  The  solution  starts  to  oscillate 
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in  time  and  it  takes  about  150  time  units  to  clearly  depart  from  the  original  DM. 
Sub-figure  3  corresponds  to  a  typical  solution  inside  the  instability  region  located  at 
ip  <  0.  Note  that  this  instability  does  not  have  a  clear  frequency  and  that  it  exhibits 
has  a  much  higher  growth  rate:  it  takes  only  40  time  units  to  completely  destroy 
the  original  DM. 

Fig  17  shows  the  stability  properties  of  family  DM2.  The  DM  in  this  family 
correspond  to  X  <->  —X  symmetric  defects,  and  there  are  stable  DM  with  a  single 
maximum  (for  Xi  <  0,  see  sub-figure  1)  and  also  with  two  maxima  (for  X\  >  0, 
see  sub-figure  3).  In  this  case  the  DM  become  unstable  as  X\  is  increased  (i.e. ,  as 
the  distance  between  the  maxima  grows)  and,  for  X\  =  2.5,  they  are  all  unstable, 
independently  of  the  value  of  ip.  If  |Ah|  is  further  increased,  then  the  defect  moves 
to  the  tail  of  the  GS  and  we  recover  again  the  stability  of  the  GS,  even  in  the  two 
maxima  case  Xi  >  0.  There  appears  to  be  again  two  unstable  regions:  one  large 
covering  all  ranges  of  ip,  and  a  small  one  near  ip  —  n/2.  In  sub- figure  2  a  typical 
solution  of  the  large  unstable  area  is  depicted;  the  growth  rate  of  the  instability  is 
quite  large  and  the  original  DM  is  completely  destroyed  for  T  =  70.  A  solution  from 
the  other  instability  region  is  represented  in  sub-figure  4.  As  it  happened  for  the 
DM1  family,  the  instability  onset  is  now  oscillatory  and  the  growth  rate  is  clearly 
much  smaller. 
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Figure  16:  Stability  of  the  DM1  family  as  a  function  of  Xi  and  0.  Red  dots  indicate  instability. 
The  solution  is  also  shown  for  three  representative  cases  (marked  with  black  dots).  Top  plots: 
spatial  profiles  of  |^4+|  (blue  line)  and  |v4“|  (red  line)  at  the  end  of  the  simulation,  and  the  profile 
of  the  DM  used  as  initial  condition  (black  line).  Bottom  plots:  temporal  evolution  of  the  norms 


36 


Figure  17:  Stability  of  the  DM2  family  as  a  function  of  X\  and  0.  Red  dots  indicate  instability. 
The  solution  is  also  shown  for  four  representative  cases  (marked  with  black  dots).  Top  plots: 
spatial  profiles  of  |7l+|  (blue  line)  and  |A“|  (red  line)  at  the  end  of  the  simulation,  and  the  profile 
of  the  DM  used  as  initial  condition  (black  line).  Bottom  plots:  temporal  evolution  of  the  norms 

II  d  I  4±  in 
II  dT  1^1,21  II- 
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Pulse  trapping  simulations 

In  this  section  we  explore  the  light  trapping  characteristics  of  two  representative 
defects  from  the  families  DM1  and  DM2  analyzed  in  the  previous  section.  We  send 
a  propagating  GS  and  simulate  numerically  its  collision  with  a  defect.  The  time 
evolution  of  the  system  (41)  is  computed  with  the  numerical  method  described  in 
the  next  section,  using  a  very  large  spatial  domain  [—100, 100]  and  with  resolution 
AX  =  .01  and  AT  =  .01. 

The  initial  condition  in  all  simulations  is  a  right  moving  GS  centered  at  A"  =  —40. 
The  moving  GS  are  a  well  known  (see,  e.g.,  [27,  23,  28,  29,  15])  family  of  propagating 
pulse  like  solutions  of  the  NLCME  (29)- (30)  that  can  be  written  in  the  form 

A+S(X,T)  =  -^H{C)e~y/2+i^)-iKlT 
AqS(X,  T)  =  VkH(Z)  ey/2+i^)-™iT 

where  £  =  k{X  —  cT )  is  a  moving  spatial  coordinate,  c  is  the  velocity  of  the  GS 

c  =  tanhy,  with  —  oo  <  y  <  +oo, 


K'y  is  its  temporal  frequency,  with 

7  =  cos 9/  coshy  =  cos^Vl  —  c2,  with  0  <  9  <  n, 


and  the  functions  H  and  cp  are  given  by 


m 


2crsinh(2y)  arctan{tan(0/2)  tanh[(coshy  sin#)/]} 
1  +  a  cosh  2 y 
sin  6 

a/1  +  a  cosh  2 y  cosh[(cosh  y  sin  6 )£  +  i0/ 2} 


+  (sinhy  cos  6*)^, 


Note  that  the  moving  GS  above  are  a  two  parameter,  c  —  7,  family  of  solutions  that, 
in  the  particular  case  of  zero  propagation  velocity  c  =  0,  reduce  to  the  standing  GS 
previously  used  to  obtain  the  DM.  In  Fig  18  we  represent  the  profile  of  a  moving 
GS  and  its  domain  of  existence:  c2  +  y2  =  1. 

The  first  set  of  simulations  corresponds  to  a  non-symmetric  defect  of  the  DM1 
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Figure  18:  Left:  moving  Gap  Soliton  at  t  =  0,  A+{A~)  thick  (thin)  line,  with  parameters  a  = 
k  =  1,  c  =  \  and  0  =  2.  Right:  domain  of  existence  of  the  GS  in  the  c  —  7  plane  (the  dot 
corresponds  to  the  GS  shown). 


family  labeled  1  in  Fig.  16.  The  parameters  in  system  (41)  for  this  case  are  k  —  1 
and  a  —  1/2,  and  the  transmission  and  reflection  coefficient  are  given  by 

ct  =  0.7270  -  0.4455i,  cr  =  -0.4516  -  0.2628i, 

which  correspond  to  the  parameters  in  eqs.  (50)-(51):  9  —  7t/4,  X\  =  1  and  ip  = 
0.5498.  The  associated  standing  nonlinear  DM  is  stable  and  is  represented  in  sub- 
figure  1  of  Fig  16.  The  temporal  frequency  of  this  DM  is  l/\/2  (see  eqs.  (42),  (43) 
and  (44))  so,  in  order  to  try  to  excite  it,  we  send  moving  GS  solitons  from  the  left 
with  the  same  temporal  frequency,  7  =  l/v^2,  and  velocities  c  =  0.2,  0.3,  0.4, 0.5,  0.6. 
The  results  are  presented  in  Figs.  19  to  23,  where  the  X  —  T  diagram  of  the  light 
intensity  during  the  collision  of  the  moving  GS  with  the  defect  (located  at  X  =  0)  is 
represented  together  with  the  profiles  of  the  solution  at  initial  and  final  time.  Note 
that  for  low  velocities  c  =  0.2  and  c  =  0.3  (Figs.  19  and  20)  the  GS  is  basically 
reflected,  and  for  higher  velocities  c  =  0.4,  c  =  0.5  and  c  =  0.6  (Figs.  21,  22  and 
23)  the  GS  just  passes  through  the  defect.  In  all  cases  some  radiation  is  shed,  but 
no  light  intensity  remains  trapped  at  the  defect.  The  speed  of  the  pulse  coming  out 
of  the  collision  with  the  defect  is  negative  for  c  =  0.3  (see  Fig.  20)  and  positive  for 
c  =  0.4  (see  Fig.  21),  so  one  can  try  to  adjust  c  to  end  up  having  a  standing  pulse 
after  the  collision.  But  this  kind  behavior  with  GS  trapping  only  for  a  single  value 
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of  c  is  a  no n- robust,  unstable  situation  that  has  not  much  interest  since  in  a  realistic 
setting  we  will  never  be  sending  the  exact  GS,  and  therefore  we  will  always  end  up 
with  a  pulse  slowly  drifting  away  from  the  defect. 

In  the  second  set  of  simulations  we  analyze  a  symmetric  defect  of  the  DM2  family. 
To  be  more  precise  the  one  that  is  labeled  1  in  Fig  17,  which  corresponds  to  k  —  1 
and  a  =  1/2,  transmission  and  reflection  coefficient 

ct  =  0.5071  +  0.7326i,  cr  =  -0.3733  -  0.2584i, 

and  parameters  6  =  7t/4,  Xi  =  —1  and  i/j  =  —0.4712  in  eqs.  (52)-(53).  The 
associated  DM  is  stable  (see  sub-figure  1  of  Fig  17)  and  has  again  temporal  frequency 
I/a/2-  The  initial  condition  in  all  simulations  is  again  a  left  going  GS,  centered  at 
X  =  —40,  with  frequency  7  =  I/a/2,  and  velocities  c  =  0.2,  0.3,  0.4,  0.5,  0.6.  The 
results  of  the  collisions  with  the  defect  are  presented  in  Figs.  24  to  28.  For  small 
velocity  c  =  0.2  (see  Fig.  24)  a  large  part  of  the  intensity  of  the  incoming  GS  is 
trapped  by  the  defect,  that  ends  up  exhibiting  a  pulsating  localized  defect  mode. 
As  the  velocity  of  the  GS  is  increased,  the  intensity  of  the  trapped  pulse  is  reduced 
and  its  oscillatory  character  ceases  at  c  =  0.4  (see  Figs.  26).  Also,  from  c  =  0.4  and 
above  an  outcoming  pulse  from  the  collision  with  the  defect  can  be  clearly  identified 
(see  Figs.  26,  27  and  28),  which  looks  like  a  GS  different  from  the  initially  sent 
one.  The  situation  is  now  for  this  symmetric  defect  completely  different:  for  all  GS 
velocities  there  is  always  trapping  of  light  intensity  by  the  defect,  in  what  appears 
to  be  a  quite  robust  behavior. 
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Figure  19:  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.2  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  \A~\  in  red  at  T  =  0  (thin  lines)  and  T  =  350  (thick  lines). 
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Figure  20:  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.3  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  300  (thick  lines). 
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Figure  21:  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.4  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  250  (thick  lines). 
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Figure  22:  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.5  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  200  (thick  lines). 
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Figure  23:  DM1  defect  interaction  with  a  GS  of  speed  c  =  0.6  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  200  (thick  lines). 
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Figure  24:  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.2  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  350  (thick  lines). 
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Figure  25:  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.3  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  300  (thick  lines). 
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Figure  26:  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.4  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  250  (thick  lines). 
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Figure  27:  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.5  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  |^4~|  in  red  at  T  =  0  (thin  lines)  and  T  =  200  (thick  lines). 
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Figure  28:  DM2  defect  interaction  with  a  GS  of  speed  c  =  0.6  sent  from  the  left.  Bottom:  X  —  T 
diagram  of  the  light  intensity,  I  =  -\/|v4+|2  +  |A_|2,  during  collision.  Top:  profiles  of  |A+|  in  blue 
and  \A~\  in  red  at  T  =  0  (thin  lines)  and  T  =  200  (thick  lines). 
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Numerical  simulation  of  light  propagation  in  a  FBG  with  a  defect 

To  simulate  the  dynamics  of  light  propagation  in  an  FBG  with  a  sharp  defect  we 
have  to  numerically  integrate  the  problem  (41)  in  a  finite  domain  [—L,L\: 


A+t  -  A+x  =  inAi  +  iAt(a\At\2  +  |Ar|2), 
A^T  +  A^X  =  i/tAf  +  %A^  (crjA^  I"  +  \Af\2), 

for 

X  <  0, 

A2t  ~  A2x  ~  ^kA2  +  iA2 (o' \A2\  +  \A2\  ), 
A2p  +  A2y  =  inA2  +  %A2  (cr  |  a!2  1  +  \At\  )j 

for 

X  >  0, 

together  with  the  jump  conditions  (39): 

A±  =  CtA2  +  crAt,  and  Af  =  ctAt  —  cr— 

■At, 

at  X  = 

ct 


and  the  appropriate  boundary  conditions  at  X  =  ±L: 

Ai(-L,T)=F1(T)  at  X  = -L  ,and  A+  =  F2(T),  at  X  —  L. 


We  use  a  uniform  spatial  grid  with  spacing  A  A"  =  L/N,  and  a  second  order 
approximation  for  the  amplitudes  and  its  spatial  derivatives  at  the  mid  points  of  the 
grid  (see  Fig  29).  The  resulting  discretized  ODEs  are  given  by 


d  ,  Af+1  +  Af , 


d  Ai+1  +  Ai 


dT 
i 

dT 


,A++1-Af  .  At  -At  NL++1+NLf 

-  AX  - -)  + 

Ar+i  -  Ai 


XX 


+  ini 


.At, -At.  nlt+1+nl; 


j+i 


+ 


+  0(  XX2), 
+  0(  XX2), 


where  the  terms  NL±  account  for  the  nonlinear  terms  of  the  equations  for  A±, 
respectively.  We  can  rewrite  the  equations  above  in  matrix  form  as 


Af 

’  1  0  1  o  " 

d 

A. ■ 

-(-H- 

Ax 

’-101  o' 

+  in 

’  1  0  1  0  ' 

0  10  1 

dT 

At  i 
Ai+l 

0  10-1 

0  10  1 

+ 


1  o 
0  1 


1  0 
0  1 


NL+ 

NL- 

NLf+1 

NLi+i 


At 
A- 
Ai+ 1 
A+ 1 
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Figure  29:  Discretization  of  the  spatial  domain  [— L,L]. 


We  apply  these  equations  at  each  mid  point  in  regions  1  and  2  (see  Fig  29)  and 
obtain  AN  equations  for  the  4 (N  +  1)  amplitudes  at  the  grid  points,  which  we  store 
in  the  vector 


A  —  [dlgi,  A1  . . . ,  A^Ar_)_1,  v41)7V+i,  A^,  A21i 


4  +  A~  C 

<  /12,7V+lJ  • 


The  remaining  4  equations  are  obtained  from  the  jump  conditions  (39), 


1  —  cr  — ct  0 

4  + 

^qjv+i 

^l,AT+l 

0 

0  -ct  cr  |  1 

0 

^2,1 

and  the  boundary  conditions 


Ai.\  —  F\  (T)  ,and  A2^n+i 


Fi(T). 


The  complete  discretized  problem  can  be  casted  as  a  system  of  4(iV+l)  differential 
algebraic  equations  of  the  form: 

dA 

Mt—  =  MlA  +  MnlNL  +  MbcA  -  BC ,  (54) 

where  the  4 (N  +  1)  x  4 (N  +  1)  matrices  MT,  ML,  MBC  and  MNL  do  not  depend  on 
T  and  are  highly  sparse,  and  the  vectors  NL  and  BC  are  4 (N  +  1)  long  and  contain 
the  nonlinear  terms, 


NL  =  [. NL+ ,,  iVLM, . . . ,  NL+n+1,  NL~n+1 ,  JVL+,,  NL^, 


<  nl2,n+1,  nl2'N+1]  , 
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and  input  from  the  boundary  conditions, 


BC=[F1(T),0,...,0,F2(T)}t. 

The  algebraic  equations  come  from  the  jump  and  the  boundary  conditions,  and  are 
placed  at  the  rows  1,  2 (TV  + 1),  2(TV  + 1)  + 1,  and  4 (TV  + 1)  of  the  system.  The  matrix 
Mbc  is  different  from  zero  only  in  these  rows,  where  the  rest  of  the  the  matrices 
Mt,  Ml,  and  MNL  do  not  contribute. 

Finally,  the  time  integration  of  (54)  is  performed  using  a  second  order  Runge- 
Kutta  scheme.  In  order  to  ensure  the  stability  of  the  resulting  scheme,  the  non¬ 
linear  terms  are  integrated  explicitly  and  the  linear  terms  (including  the  algebraic 
equations)  are  implicitly  treated  (see  [30]). 

The  first  sub-step  gives  an  estimate  A*  of  the  solution  at  T  +  AT : 

AT 

Mt(A *  -  AT )  =  —Ml(A*  +  AT )  +  A TMnlNLt  +  MBCA*  -  BC(T  +  AT), 

and  with  this  estimate  we  correct  the  integration  of  the  nonlinear  term  in  the  second 
and  final  sub-step: 

Mt(At+at  -  At )  =  ^-Ml{At+at  +  At )  +  ^Lmnl(NLt+at  +  NLt ) 

+  MbcAt+at  -BC(T  +  AT). 

Note  that  at  every  sub-step  we  have  to  solve  the  same  linear  system.  This  allows 
us  to  factorize  the  system  matrix  once  at  the  beginning  of  the  integration  and  then, 
at  each  sub-step,  we  have  only  to  solve  the  factorized  system.  The  simulation  code 
has  been  programmed  in  MATLAB,  and  runs  quite  fast;  it  takes  about  30  seconds 
in  a  standard  3GHz  Pentium  PC  to  advance  10000  time  steps  in  a  grid  with  2000 
spatial  points. 

In  order  to  check  that  the  code  is  actually  second  order  in  time  and  space  we  run 
the  following  test:  in  a  domain  with  L  =  4,  and  with  ct  —  \  +  |  and  cr  —  \  —  |,  we 
send  a  pulse  from  the  left  side, 

Fl(T)  =  e-(T~4)2  and  F2(T)  =  0, 

and  integrate  the  system  with  zero  initial  condition  from  T  =  0  to  T  =  10  (the 
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final  solution  at  T  =  10  is  quite  involved  as  it  can  be  seen  in  Fig  29).  We  perform 
the  integration  for  three  different  temporal  and  spatial  resolutions  AX  =  AT  = 
.01,  .01/2  and  .01/4.  In  Fig  31  we  plot,  using  logarithmic  scale,  the  differences 
between  two  successive  solutions  versus  the  spatial  coordinate  X  at  T  =  10.  It  can 
be  clearly  appreciated  that  the  ratio  between  the  differences  is  4,  which  indicates 
that  the  solution  is  second  order,  i.e. ,  the  error  goes  as  AX2  +  AT2. 


Figure  30:  Spatial  profiles  of  the  solution  at  T=10:  |>1+|  (blue  line)  and  |d+|  (red  line). 


Figure  31:  Modulus  of  the  difference  between  two  increasingly  refined  solutions  at  T=10.  Blue  line 
for  A+and  red  for  A~ . 
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Conclusions 

In  this  report  we  have  analyzed  the  weakly  nonlinear  dynamics  of  nonlinear  light 
propagation  in  a  FBG  with  a  narrow  defect.  The  problem  has  been  first  formulated 
in  terms  of  two  sets  of  NLCME  that  apply  at  both  sides  of  the  defect  (regions  1  and 
2  in  Fig.  6)  together  with  a  jump  condition  at  the  location  of  the  defect,  see  eqs. 
(41).  The  coefficients  for  the  jump  condition  depend  on  the  profile  of  the  defect,  and 
their  computation  requires  to  solve  a  scattering  problem.  This  problem  has  been 
solved  numerically  for  several  representative  defects  in  order  to  illustrate  the  relation 
between  the  defect  shape  and  the  resulting  coefficients.  From  the  standing  GS  exact 
solution  we  have  constructed  a  general  family  of  nonlinear  defect  modes  that  depend 
on  the  particular  values  of  the  coefficients  of  the  jump  condition.  The  defect  modes 
correspond  to  persistent  in  time  light  distributions  that  remain  spatially  confined 
to  the  vicinity  of  the  defect.  The  stability  of  these  nonlinear  defect  modes  has  been 
analyzed  numerically,  and  for  two  representative  stable  defect  modes  we  have  also 
explored  the  possibility  of  light  pulse  trapping  (i.e.,  the  possibility  of  exciting  the 
standing  defect  modes  by  sending  moving  pulses  to  collide  with  them).  A  numerical 
code  for  the  simulation  of  the  temporal  evolution  of  the  complete  problem  of  light 
propagation  in  a  FBG  with  a  narrow  defect,  eqs. (41),  has  been  also  implemented, 
and  its  details  are  described  in  this  report.  This  numerical  code  is  required  for  the 
study  of  the  defect  mode  stability,  and  for  the  determination  of  the  light  trapping 
characteristics  of  the  defects. 

Some  conclusions  that  can  be  drawn  from  this  research  project  are  now  summa¬ 
rized  below: 

1.  A  narrow  defect  gives  rise  to  a  first  order  effect  only  if  the  distortion  of  the 
refraction  index  that  it  produces  is  comparable  to  the  refraction  index  itself. 
In  terms  of  the  scaling  of  MLE  eqs.  (15)  and  (16),  the  size  of  the  defect  has 
to  be  of  order  1,  much  larger  that  the  amplitude  of  the  grating,  which  is  of  the 
order  of  An  <C  1. 

2.  If  the  narrow  defect  has  small  amplitude,  of  the  order  of  the  grating  amplitude, 
then  its  effect,  in  first  approximation,  is  simply  not  felt.  This  looks  surprising 
at  first,  but  it  makes  perfect  sense  because  if  it  were  not  so,  every  imperfection 
present  in  the  small  grating  would  produce  a  first  order  effect  and,  since  there 
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are  always  imperfections  clue  to  manufacturing  and  wear,  it  would  be  impossible 
in  practice  to  see  any  of  the  dynamics  associated  with  periodic  gratings. 

3.  The  effect  of  a  narrow  defect  is  a  jump  of  the  envelopes  along  the  characteristic 
lines  of  the  hyperbolic  NLCME.  This  has  been  previously  modeled  by  adding 
Dirac-delta  functions  to  the  NLCME  [16,  17,  18].  This  formulation  covers 
only  the  particular  case  of  x-symmetric  defects  and,  even  in  this  case,  it  does 
not  provide  information  about  the  relation  between  the  deffect  shape  and  the 
resulting  jump. 

4.  A  family  of  nonlinear  defect  modes  has  been  obtained  by  joining  pieces  of 
standing  GS.  This  procedure  was  already  used  in  [17],  but  the  family  that  they 
derived  was  much  more  reduced  because  it  only  covered  the  particular  case 
of  symmetric  defects  with  real  transmission  coefficient.  We  also  show  in  this 
report  that  these  more  general  defect  modes  that  were  not  considered  in  [17] 
can  be  stable. 

5.  The  formulation  presented  in  this  paper  in  terms  of  two  sets  of  NLCME  plus 
some  jump  conditions  (see  eqs.(41))  is  also  better  suited  for  numerical  com¬ 
putations  than  the  previous  Dirac-delta  formulation  used  in  [16,  17,  18].  The 
Dirac-deltas  have  to  be  approximated  by  introducing  very  high  values  of  the 
functions  at  a  very  narrow  region  of  the  spatial  grid,  and  this  procedure  is  ob¬ 
viously  problematic  from  the  point  of  view  of  ensuring  the  order  of  accuracy  of 
the  numerical  method.  Our  formulation  has  no  problems  of  this  kind  since  it  is 
just  composed  of  two  sets  of  perfectly  regular  PDE’s  (which  can  be  discretized 
in  a  completely  straight  forward  way)  together  with  some  coupling  boundary 
conditions  that  account  for  the  jumps  of  the  amplitudes  at  the  defect. 

6.  The  light  pulse  trapping  characteristics  of  the  defects  have  been  also  explored 
and,  generally  speaking,  it  has  been  found  that  while  non-symmetric  defects 
appear  to  just  deflect  the  pulses  without  trapping  any  light,  the  symmetric 
defects  always  end  up  producing  states  with  persistent  trapped  light  at  the 
defect. 

Finally,  it  is  interesting  to  mention  that  the  formulation  and  the  tools  developed  in 
this  research  project  will  surely  be  very  useful  for  future  studies  on  this  problem.  In 
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particular,  they  can  be  used  to  address  a  question  that  is  very  important  from  the 
technological  point  of  view:  how  to  design  the  internal  profile  of  the  defect  in  order 
to  produce  some  prescribed  light  reflection  or  trapping  properties. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


x 

An 

AT, AX 

Sd 

7 

K 

\ a 
Vo 

uk,  U) 

Qp,  UJp 
<J 

£o 

a,  a',  ai,a2 
A±,Af,At 

B 


c 

Ct,  Cr,  Cf  ,  cr 

D 

E 

fix) 

k 

L 


no,  n 
P 


x,  t,  X,  T 


Xoo 


DM 

DM1,  DM2 
FBG 


linear  polarizability 

grating  strength 

time  step,  spatial  grid  spacing 

length  of  the  defect 

cubic  Kerr  effect  coefficient 

rescaled  grating  strength  in  the  NLCME 

grating  period 

vacuum  permeability 

propagating  light  frequency 

characteristic  frequency  of  the  polarization  response 
rescaled  nonlinear  coefficient  in  the  NLCME 
vacuum  permittivity 

solutions  of  the  defect  scattering  problem 

slow  envelopes  of  the  light  wavetrains 

magnetic  field 

vacuum  speed  of  light 

transmission  and  reflection  coefficients 

dielectric  displacement 

electric  field 

defect  profile 

propagating  light  wavenumber 

spatial  domain  [— L,  L\  for  the  numerical  integration 

bare  fiber  and  FBG  refractive  index 

polarization  field 

propagating  light  group  velocity 

spatial  and  temporal  variable,  fast  and  slow 

finite  length  for  the  numerical  integrations 

Defect  Mode 

Defect  mode  families 

Fiber  Bragg  Grating 
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GS 

MLE 

NLCME 


Gap  Soliton 

Maxwell  Lorentz  Equations 

Non  Linear  Coupled  Mode  Equations 
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